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ABSTRACT 

A new numerical code, designed for the detailed numerical treatment of nonlinear diffusive 
shock acceleration, is used for modeling of particle acceleration and radiation in young supernova 
remnants. The model is based on spherically symmetric hydrodynamic equations complemented 
with transport equations for relativistic particles. For the first time, the acceleration of elec- 
trons and protons by both forward and reverse shocks is studied through detailed numerical 
calculations. We model the energy spectra and spatial distributions of nonthermal emission of 
the young supernova remnant RX J1713. 7-3946 and compare the calculations with the spectral 
and morphological properties of this object obtained in broad energy band from radio to very 
high energy gamma-rays. We discuss the advantages and shortcomings of the so-called hadronic 
and leptonic models which assume that the observed TeV gamma-ray emission is produced by 
accelerated protons and electrons, respectively. We discuss also a "composite" scenario when 
the gamma-ray flux from the main parts of the shell has inverse Compton origin, but with a 
non-negligible contribution of hadronic origin from dense clouds interacting with the shell. 



Subject headings: cosmic rays- acceleration- instabilities 



1. Introduction 

The paradigm of the diffusive shock accelera- 
tion (DSA) of relativ istic particl es (Krvmsky ll977 



Axfordet al. Il977l Bell Il978l Blandford & Os 



triker l"l978h is generally accepted as the most likely 
scenario of production of galactic cosmic rays (CR) 
in supernova remnants (SNRs). Over the last 30 
years a significant progress has been achieved in 
the development of theoretical models and under- 
standing of the basic features of DSA (see e.g. 
Malkov & Drurv l2001 for a review). On the other 
hand, the recent detailed studies of spectral and 
morphological features of young SNRs, first of all 
in the X-ray and very-high-energy or TeV gamma- 
ray band, provide excellent observational mate- 



rial for development of detailed numerical models 
of acceleration and radiation of relativistic elec- 
trons and protons in young SNRs. These observa- 
tions generally confirm in general terms the pre- 
dictions of DSA. In particular the synchrotron X- 
radiation observed from many young SNRs im- 
plies an existence of multi-TeV energy electrons 
which is naturally explained by DSA. The detec- 
tion of TeV gamma-rays from SNRs, like Cas A, 
RX J1713 7-394 6, Vela Jr., RCW 86 (see Aharo- 
nian et al. 120081 for a recent review) give a more di- 
rect and unambiguous information about the effec- 
tive acceleration of particles, electrons and/or pro- 
tons, in SNRs to energies exceeding 100 TeV. In 
spite of the high quality data obtained from these 
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SNRs and intensive theoretical and phenomeno- 
logical studies, we cannot yet firmly distinguish 
between two competing processes of gamma-ray 
production: the inverse Compton (1C) scattering 
of the background electromagnetic radiation by 
accelerated electrons (leptonic model) or a decay 
of neutral pions, mainly produced by accelerat- 
ing protons during the interaction with the gas 
of the remnant (hadronic model). It should be 
noted that this ambiguity concerns the interpre- 
tation of the origin of gamma-rays, rather than 
the question of acceleration of protons in general. 
The presence of multi-TeV electrons in the SNR 
shells most likely implies also acceleration of pro- 
tons - not only because the DSA operates iden- 
tically for both electrons and protons; the pres- 
ence of larger amount of nonthermal protons is 
needed in any case for production of magnetohy- 
drodynamic (MHD) turbulence in the vicinity of 
a supernova shock. This is a necessary condition 
for effective acce leratio n of particles via DSA (see 
Malkov & Drury l200ll for a review) . 

In this paper we conduct detailed study accel- 
eration of electrons and protons with an empha- 
sis on the spectral and morphological features of 
high energy radiation produced by these particles 
in young supernova remnants. For that purpose 
we use a new numerical code of nonlinear dif- 
fusive shock acceleration developed by one of us 
in colla borati on with V. Ptuskin (Zirakashvili & 
Ptuskin[200i). This model can be considered as 
a natural development of exist ing numerica l code s 
(see e.g. Bcrezhko et al ll994 Kang et al 120061 1 . 
with new additional elements which despite their 
strong impact on the overall picture of acceleration 
in general, and on the properties of high energy ra- 
diation of SNRs in particular, have been ignored in 
the past. Namely, in our treatment, the solution 
of spherically symmetric hydrodynamic equations 
is combined with the energetic particle transport 
and acceleration by the forward shock (FS) and 
reverse shock (RS). Nonlinear response of ener- 
getic particles via their pressure gradient results 
in self-regulation of acceleration efficiency To our 
knowledge, this is the first study in which the mod- 
ification of the reverse shock and the nonthermal 
radiation related to the reverse shock is taken into 
account. The detailed calculations of radio, X- 
ray, gamma-ray emission components conducted 
within a self-consistent treatment of particle ac- 



celeration by both forward and reverse shocks al- 
lows direct comparison of the observed spectral 
and morphological features with the model predic- 
tions. The inclusion of the radiation components 
related to the reverse shock seems to us a rather 
obligatory condition, given the effective accelera- 
tion of particles at the reverse shock to very high 
energies, as it is demonstrated in this paper. In 
this regard we note that the parameters charac- 
terizing the reverse shock can be significantly dif- 
ferent in comparison with the parameters at the 
forward shock. As a result, the properties of ra- 
diation components from the reverse and forward 
shocks can be also significantly different. In par- 
ticular the density of plasma and the magnetic 
field in the reverse shock can be very low which 
would dramatically increase the contribution of 
the IC component compared to the hadronic (tt°- 
decay) component of gamma-rays. Because of the 
higher gas density and stronger magnetic field, in 
the forward shock one may expect an opposite re- 
lation between the contributions of the electrons 
and protons to the high energy gamma-radiation. 

The results of this study have a general char- 
acter and can be applied to different young SNRs. 
In order to demonstrate their astrophysical impli- 
cation, in this paper we apply the model to the su- 
pernova remnant RX J1713. 7-3946 - the brightest 
and best studied TeV gamma-ray source among 
known shell-t ype supernova remnants (Aharonian 



et al. l2007af ). The high quality gamma-ray and 



X-ray images and energy spectra available for this 
source make it a kind of template source for the- 
oretical studies. On the other hand, this source 
have quite unusual radiation features (e.g. a lack 
of thermal X-radiation, a very low intensity of the 
synchrotron radio emission, etc.) so the conclu- 
sions derived for this source cannot be directly ex- 
tended to other SNRs. 

The paper is organized as follows. In Section 
2 we briefly summarize the observational informa- 
tion about the supernova remnant RX J1713.7- 
3946. The short description of the model and basic 
features of accelerated particles is given in Section 

3. The results of modeling of the broad-band emis- 
sion radiation produced in hadronic and leptonic 
scenarios are presented and discussed in Sections 

4, 5 and 6. In Section 7 we compare the magnetic 
field amplification and the maximum particle en- 
ergies achieved in forward and reverse shocks. Fi- 
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nally, in Sections 8 and 9 we discuss and summa- 
rize the obtained results. 

2. Observational properties of RX J1713.7- 
3946 

The shell-type SNR RX J1713. 7-3946 was 
discovered in X-rays during the ROSAT All- 
Sky Survey (Pfeffermann & Aschenbach 19961 ). 
The following observations of ASCA, Cha ndra , 
XMM-Newton and Su zaku (Koyama et al. 1997 , 
Uchiyama et al. 20031 Cassam-Chen ai et a l. .2004, 
Hiraga et al. 120051 Takahashi et al. 120081 Tanaka 
et al. 2008) revealed that the X-ray emission 
predominantly consists of the synchrotron com- 
ponent without an indication of a thermal X-ray 
component at the level of sensitivity of these in- 
struments. The Suzaku measurements (Tanaka 



2001 



ct al. 20081 ) allowed accurate derivation of the 



broad-band X-ray spectrum of the source over 
two energy decades, from 0.4 keV to 40 keV. It 
is described by a power-law with a photon index 
~ 2 and a smooth cut-off around a few keV. 

The remnant is close to the position of the 
guest star AD393 de scribe d in the ancient Chinese 
records (Wang et al ll997l h If this not an acciden- 
tal coincidence, the age of the remnant should be 
close to 1600 years. 

Using the association of the supernova remnant 
with nearby molecular clouds, different estimates 
of the distance to the so urce h ave been suggested: 
D = 6 kpc (Slane et al. Il999h . D = 1 kpc (Fukui 
et al. 2003f ) and 1.3 ± 0.4 kpc (Cassam-Chenai et 
al. l2004h . 

The remnant was claimed to be detected in 
gamma-rays by th e CA NGAROO collaboration 
(Muraishi et al. l2000h . Later, the observa- 



tion with the HESS array of imaging atmospheric 
Cherenkov telescopes revealed indeed a bright TeV 
gamma-ray source with a shell-type morphology 
quite similar to the synchrotr on X- ray image of 



the source (Aharonian et al. 120041 ) The energy 
spectrum measured over the interval from 0.3 TeV 
to 100 TeV is hard; it is characterized by a photon 
index ~ 2 with a break o r a cutoff around 10 TeV 
(Aharonian et al. l2007al ). 

The remnant is rather faint in radio. The re- 
cent ATCA measurements at 1.4 GHz gave the 
radio-intensity 22 ± 2 Jy from the whole remnant 
(Acero et al. 20091 ). The radio and X-ray images 



clearly show als o an in ner shell (Ellison et al 
Lazendic et al. 2004 Cassam-Che nai et al. 
Hiraga et al. I2005L Acero et al. l2009h with an 
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angular diameter 0.5 degrees while the angular di- 
ameter of the remnant is close to one degree. The 
inner shell in the radio and X-ray images can be 
naturally attributed to the reverse shock propa- 
gating in the supernova ejecta, assuming that the 
electrons are accelerated also at the reverse shock. 
This is likely to be the case of the supernova rem- 
nant Cas A (Helder & Vink l2008l ). Generally, the 
reverse shock is not treated as an efficient accel- 
erator, because the magnetic field of ejecta might 
be very weak due to the large expansion factor of 
the exploded star. However, similar to the case of 
the forward shock, the magnetic field can be sig- 
nificantly amplified also at the reverse shock (see 



Ellison et al. 120051 ). for example, in the course of 



the no nreson ant streaming instability introduced 
by Bell [2004|. 

Another probable explanation of the inner shell 
could be its association with a pulsar wind nebula. 
It is known that an active pulsar may produce a 
corresponding ne bula in an expanding supernova 
ejecta (Chevalier [2005). The inner shell does con- 
tain the radio pulsar PSR J1713-3945 with a pe- 
riod P = 392 ms, but the latter is located at a 
distance D = 4.3 kpc. Also, the pulsar's cur- 
rent spin-down luminosity does not exceed 10 
erg s _1 . This power is too low to produce the in- 
ner shell. The shell contains also a central compact 
object 1WGA J1713.4-3949 which is presumably 
a neu tron star related to the re mnant (Lazendic 
ct al. [20031 Cassam-Chenai et al. I2004J ) 



3. Nonlinear model of diffusive shock ac- 
celeration 

Hydrodynamical equations for the gas density 
p(r,t), gas velocity u(r, i), gas pressure P g (r,t), 
and the equation for the quasi-isotropic CR mo- 
mentum distribution N(r,t,p) in the spherically 
symmetrical case are given by 



dp 



1 d 



— — = T Tp r U P 



at 



r 



dr 



(1) 



— - - u — - - (^-2- + — ) r'l 

dt dr p \ dr dr ' 



3 



dt 

ON 
~dt 



dPg 

dr 



jP g dr 2 v 
r 2 dr 



( % -l)(w-u) — 
(3) 



_L_5 

r 2 <9r 



r 2 D(p, r, t) 



dN dN dN p dr 2 w 



dr 



dr dp 3r 2 dr 



4-7rp 2 m 

r/ b 5(p-p b ) 
4np 2 m 



p(Rf+0, t){R f -u{R+0, t))S(r-R f (t)) 



p(Rb-0, t)(u(R b -0, t)-R b )6{r~R b {t)) 



(4) 

Here P c = Air f p 2 dpvpN/3 is the CR pressure, 
w(r, t) is the advective velocity of CRs, j g is the 
adiabatic index of the gas and D(r, t,p) is the CR 
diffusion coefficient. It was assumed that the diffu- 
sive streaming of CRs results in the generation of 
magnetohydrodynamic (MHD) waves. CR parti- 
cles are scattered by these waves. That is why the 
CR advective velocity w may differ from the gas 
velocity u. Damping of these waves results in an 
additional gas heating. It is described by the last 
term in Eq. (3). Two last terms in Eq. (4) cor- 
respond to the injection of thermal protons with 
momenta p — pf, p — p b and mass m at the fronts 
of the forward and reverse shocks at r = Rf(t) 
and r = R b (t) respectively. The dimensionless pa- 
rameters r\f and m, determine the corresponding 
injection efficiencjtJ. 

The pressure of energetic electrons is neglected, 
i.e. the electrons are treated as test particles. The 
evolution of electron distribution is described by 
an equation similar to Eq. (4) , but with additional 
terms describing synchrotron and IC losses. 

CR diffusion is determined by magnetic inho- 
mogeneities. Strong streaming of accelerated par- 
ticles changes medium properties in the shock 
vicinity. In particular, the CR streaming insta- 
bility in young SNRs resul ts in a high level of 
MHD turbulence (Bell ll978j) and even in amplifi- 
cation of magnetic field (Bell l2004 ). Due to this ef- 
fect the protons can be accelerated to energies be- 
yond the the so-called Lagage and Cesar sky limit, 
E < 100 TeV (Lagage & Cesarsky [l983|). 



1 Throughout the paper we use the subscripts "f and "b" 
to the parameters characterizing the forward and reverse 
(backward) shocks, respectively. 



According to the recent numerical modeling of 
this instability, the magnetic field is amplified by 
the flux of run-away highest energy particles in 
the relatively broad regi on up stream of the shock 
(Zirakashvili & Ptuskin 120081 ). Magnetic energy 
density is a small fraction (~ 10~ 3 ) of the energy 
density of accelerated particles. This amplified al- 
most isotropic magnetic field can be considered as 
a large-scale magnetic field for lower energy parti- 
cles which are concentrated in the narrow region 
upstream of the shock. Streaming instability of 
these particles produces MHD waves propagating 
in the direction opposite to the CR gradient. This 
gradient is negative upstream of the forward shock 
and MHD waves propagate in the positive direc- 
tion. 

We apply a finite-difference method to solve nu- 
merically Eqs. (1-4) upstream and downstream of 
the forward and reverse shocks. The auto-model 
variables £i = r/Rf(t) and £4 = r/R b (t) are used 
instead of the radius r upstream of the forward 
shock at r > Rf and upstream of the reverse shock 
at r < R b respectively. The gases compressed at 
forward and reverse shocks are separated by a con- 
tact discontinuity (CD) at r = R c that is situ- 
ated between the shocks. We use variables £2 = 
(r - R c )/(R f - R c ) and £3 = (r - R C )/(R C - R b ) 
instead of r between the shocks downstream of the 
forward and reverse shocks, respectively. 

A non-uniform numerical grid upstream of the 
shocks at r > Rf and r < R b allows to resolve 
small scales of hydrodynamical quantities appear- 
ing due to the pressure gradient of low-energy 
CRs. Eq. (4) for CRs was solved using an implicit 
finite-difference scheme. An ex plicit conservative 
TVD scheme (Trac & Pen l2003l) for hydrodynam- 
ical equations (1-3) and uniform spatial grid were 
used between the shocks. These equations are 
solved in the upstream regions using an implicit 
finite-difference scheme. 

The magnetic field plays no dynamical role in 
the model. Since we do not model the amplifica- 
tion and transport of magnetic field here, its co- 
ordinate dependence should be specified for deter- 
mination of cosmic ray diffusion and for calcula- 
tion of synchrotron emission and losses. We shall 
assume below that the coordinate dependencies of 
the magnetic field and the gas density coincide up- 
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stream and downstream of the forward shock: 



B(r) = V 47r Po 



RfP 

M f A p 



r > R c . 



(5) 



Here po is the gas density of the circumstellar 
medium. The parameter M A is similar to the 
Alfven Much number of the shock and determines 
the value of the amplified magnetic field strength. 
The magnetic energy downstream of the shock is 
estimated to be several percents of the dynami- 
cal pressure poR 2 as was derived from the width 
of X -ray filaments in young SNRs (Volk et al. 
2005). The characteristic range of this parame- 



ter is M J A ~ 10 -j- 40. For example, the energy 
of the magnetic field of about 3. 5 per cent of the 



dynamical pressure (Volk et al. 120051 ) and char- 



acteristic compression ratio of a modified shock 
a = 6 correspond to M A ss 23. Since the plasma 
density p decreases towards the contact disconti- 
nuity downstream of the forward shock, the same 
is true for the magnetic field strength according 
to Eq. (5). This seems reasonable because of a 
possible magnetic dissipation in this region. 

Situation is different downstream of the reverse 
shock at R\, < r < R c . The plasma flow is as a 
rule strongly influenced by the Rayleigh- Taylor in- 
stability that occurs in the vicinity of the contact 
discontinuity and results in generation of MHD 
turbulence in this region. We assume that the 
magnetic field does not depend on radius down- 
stream of the reverse shock while the dependence 
in the upstream region is described by the equa- 
tion similar to Eq. (5): 

-\R b - u(r m )\ 



B(r) = y^pl 



1, r < r m , 

p/pm, r m <r < R b , 

p(R b + 0)/p m , Rb<r <R C 



(6) 



Here r m < Rb is the radius where the ejecta den- 
sity has a minimum and equals p m . This radius 
r m is generally close to the reverse shock radius 
Rb and is equal to it if the reverse shock is not 
modified by the cosmic ray pressure. 

CR advective velocity may differ from the gas 
velocity on the value of the radial component of 
the Alfven velocity Vat calculated in the isotropic 
random magnetic field: w = u + S,aVat- Here 



the factor £4 describes the possible deviation of 
the cosmic ray drift velocity from the gas velocity. 
Using Eq. (5) we obtain 



M{ V 3/3 



(7) 



The similar expression for the cosmic ray drift ve- 
locity is used upstream of the reverse shock at 
r < Rb- We use values £4 = 1 and £a = — 1 
upstream of the forward and reverse shocks re- 
spectively, where Alfven waves are generated by 
the cosmic ray streaming instability and propa- 
gate in the corresponding directions. The damp- 
ing of these waves heats the gas u pstream of the 
shocks (see McKenzie & Voik ll982l ) and limits the 
total compression ratios by a number close to 6. 
In the downstream region of the forward and re- 
verse shock at Rb < r < Rf we put £a — and 
therefore w — u. 

Here we use the energy dependence of the CR 
diffusion coefficient like in the case of Bohm diffu- 
sion: D = i]bDb calculated for the magnetic field 
radial dependencies given by Eq. (5) and (6). The 
parameter r\B describes the possible deviations of 
diffusion coefficient from the one achieved in the 
regime of Bohm diffusion, Db = vpc/3qB. Since 
the highest energy particles are scattered by small- 
scale magnetic fields, their diffusion is fast er tha n 
the Bohm diffusion (Zirakashvili & Ptuskin Exjj ). 
The same is true for low energy particles because 
they can be resonantly scattered only by a fraction 
of the magnetic spectrum. Throughout the paper 
we use the value t\b = 2. 

In real situations the level of MHD turbulence 
drops with distance upstream of the shock, so the 
diffusion could be quite fast there. The character- 
istic diffusive scale of highest energy particles is 
a small fraction £0 << 1 of t he shock radius (see 
Zirakashvili & Ptuskin l2008h and is determined 
by the generation and transport of MHD turbu- 
lence i n the upstream region ( see als o Vladimirov 
et al. 120061 . Amato & Blasi l2006h . The value 
£0 ~ In -1 (D c /D s ) is determined by ratio of dif- 
fusion coefficient D c in the circumstellar medium 
and diffusion coefficient D s << D c in the vicin- 
ity of the shock. The MHD turbulence is ampli- 
fied exponentially in time before the shock arrival 
from the background level by cosmic ray stream- 
ing instability. Since this process is not modelled 
here, we simply multiply the CR diffusion coeffi- 
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cient D to the additional factor exp((£i — l)/£o) 
upstream of the forward shock and to a similar fac- 
tor exp((l — £4)/£o) upstream of the reverse shock. 
The characteristic r ange o f £o is 0.05 -j- 0.1 (Zi- 
rakashvili & Ptuskin [2008h . 

It is believed that the supernova ejecta has some 
velocity distribution P(V ) just after the supernova 
explosion (e.g. Chevalier 1 1982f ) 



P(V) 



3(fc- 3) 



l,V< V ej 

(v/v ej y k , 



V > v f 



(8) 



>'j ■ 



Here the index k characterizes the steep power-low 
part of this distribution. The radial distribution of 
ejecta density is described by the same expression 
with V = r/t. The characteristic ejecta velocity 
V e j can be expressed in terms of energy of super- 
nova explosion E$n and ejecta mass M e j as 



f lO{k~h)E S N 



\ 3(fc - 3)M, 



1/2 



(9) 



Figures 1-5 illustrate the numerical results that 
are obtained for the SNR shock propagating in the 
medium with hydrogen number density = 0.09 
cm' 3 and temperature T — 10 4 K. The number 
density of helium nuclei O.lnn was assumed. We 
use the ejecta mass M e j = 1.5M Q , the energy of 
explosion E$n = 2.7 • 10 51 erg and the parame- 
ter of ejecta velocity distribution k = 7. These 
combination of these parameters is chosen to ex- 
plain the obtained gamma-ray fluxes within the 
hadronic model of gamma-rays (see next Section). 

The initial forward shock velocity is Vq — 4.6 • 
10 4 km s _1 . The injection efficiency is taken to 
be independent of time rjb = rjt = 0.01, and the 
injection momenta are pj = 2m(Rf — u(R + 0, t)), 
Pb = 2m(u(Rt ) — 0, fy — Rf,). The high injection effi- 
ciency results in significant shock modification al- 
ready at early epochs of the SNR expansion while 
the thermal sub-shock compression ratio is close to 
2.5. This is in agreement wit h calcu lations of col- 
lisionless shocks (Zirakashvili |2007| ) and are sup- 
ported by radio-obse rvatio ns of young extragalac- 
tic SNRs (Chevalier [2006). In hadronic models 
of gamma-radiation of RX J1713. 7-3946 the ra- 
tio of relativistic electrons to protons K ep appears 
significantly lower compared to the observed ra- 
tio for galactic cosmic rays which is around 0.01 
at 10 GeV. A possible explanation could be the 
dependence of the electron injection efficiency on 
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Fig. 1. — Time dependencies of parameters char- 
acterizing the forward and reverses shocks: the 
forward shock radius Rf (thick solid line), the 
reverse shock radius Rb (thick dashed line), the 
forward shock velocity Rf (thin solid line); the 
reverse shock velocity Rb (thin dashed line); the 
magnetic field strength downstream of the forward 
shock (thin dotted line); the ratio of the CR en- 
ergy to the total energy of the supernova explosion 
E cr /EsN (dotted line). 



the shock speed. In the paper we assume that the 
electron injection is inversely proportional to the 
square of the shock speed. Then the electron to 
proton ratio K ep ~ 10~ 4 for a characteristic shock 
speed 3000 km s" 1 of young SNRs corresponds to 
K ep ~ 10~ 2 in the old SNRs with a characteristic 
shock speed 300 km s _1 . The old SNRs probably 
produce the main part of galactic GeV electrons. 

The time evolution of the shock radii Rf and 
Rb, the forward and reverse shock speeds Vf = Rf 
and Vb = Rb, the CR energy E cr /EsN and the 
magnetic field strength B f downstream of the for- 
ward shock are shown in Fig.l. The calculations 
were performed until the present epoch at t — 1620 
yr, when the forward shock speed drops down to 
Rf — 2.76 • 10 3 km s' 1 and the forward shock ra- 
dius Rf — 10.5 pc corresponding to the angular ra- 
dius 30' of the remnant RX J1713. 7-3946, assum- 
ing 1.2 kpc distance to the source. At early epochs 
of SNR evolution the distance between reverse and 
forward shocks is only 10% of the remnant radius. 
Note that the automodel Chevalier-Nadezhin so- 
luti on wi th k — 7 predicts 23% thickness (Cheva- 
lier 1982). This can be attributed to the strong 
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Fig. 2. — Radial dependencies of the gas density 
(thick solid line), the gas velocity (dotted line), 
CR pressure (thick dashed line), the gas pressure 
(dashed line) and the electron temperature (thin 
solid line) at t = 1620 yr. The calculations result 
in the following parameters in the present epoch: 
the forward shock velocity 2760 km s _1 , its radius 
10.5 pc, the magnetic field strength downstream 
of the forward shock 127 p:G. In the same figure 
we show the positions of the forward and reverse 
shocks, (FS and RS, respectively) and the contact 
discontinuity (CD). 



modification of both shocks by CR pressure. The 
reverse shock is strongly decelerated only when 
the forward shock sweeps the gas mass compara- 
ble to the ejecta mass at t > 100 yr and when 
the transition to the Sedov phase begins. It is 
important that the observable ratio Rb/Rf ~ 0.5 
is achieved at the Sedov phase when the forward 
shock speed is close to 1/3 of the characteristic 
ejecta speed V e j. Since the synchrotron X-ray 
emission of SNR RX J1713. 7-3946 extends well 
beyond 1 keV, the present forward shock speed 
cannot be significantly less than 3 • 10 3 km s _1 . 
Therefore the characteristic velocity of the ejecta 
cannot be smaller than 10 4 km s _1 . Such high ve- 
locities may be attributed to Ib/c and lib core col- 
lapse supernovae with low ejecta masses, but not 
to the most frequent core collapse IIP supernovae 
with characteristic ejecta velocities 3 • 10 3 — 4 • 10 3 
km s^ 1 and high ejecta masses. 

Radial dependencies of several key parameters 
at the present epoch t — 1620 yr are shown in 



0.1 



0.01 



0.001 



\\\ ^N(p)p 4 c/p R 2 


P 


Vx>c: — ^nr/V 


\ \ e / V i 
















p/mc \ 
1 1 1 





0.01 



100 



10 4 



10 b 



Fig. 3. — The energy distributions of accelerated 
protons (thick lines) and electrons (multiplied to 
the factor of 5000, thin lines) at the epoch t = 1620 
yr. The spectra at both the forward shock (solid 
lines) and at the reverse shock (dashed lines) are 
shown. 



Fig. 2. In the same figure we show the positions of 
the contact discontinuity and the forward and re- 
verse shocks. The contact discontinuity between 
the ejecta and the interstellar gas is located at 
r = R c = 7.8 pc. The reverse shock in the ejecta 
is situated at r = Rb = 5.5 pc. At the Sedov stage 
the reverse shock moves in the negative direction 
and reaches the center three thousand years af- 
ter the supernova explosion. It should be noted 
that our one-dimensional calculations cannot ade- 
quately describe the development of the Raylcigh- 
Taylor instability of the contact discontinuity. In 
real situations the supernova ejecta and the cir- 
cumstellar gas are mixed by turbulent motions in 
this region (see e.g. MHD modeling of Jun & Nor- 



1996). 



In SNRs the plasma is heated to keV tempera- 
tures. The thermal X-ray emission is determined 
by the electron temperature. In young SNRs ther- 
mal electrons are not in equilibrium with pro- 
tons if the temperature equilibration proceeds 
through the minimum available electron heating, 
i.e. through Coulomb collisions (Spitzer Il968h . 
The evolution of electron temperature T e is de- 
scribed by the following equation 
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Fig. 4. — The energy distributions of accelerated 
protons (thick lines) and electrons (multiplied to 
the factor of 5000, thin lines) at t = 100 yr. Spec- 
tra at both the forward shock (solid lines) and the 
reverse shock (dashed lines) are shown. 



Fig. 5. — Spatially integrated spectra of acceler- 
ated protons (solid line) and electrons (multiplied 
to the factor of 5000, thin solid line) at t = 1620 
yr. Spectrum of run-away particles which have left 
the remnant is also shown (dashed line). 
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where X ep ~ 30 is the Coulomb logarithm. The 
corresponding temperature is also shown in Fig. 2. 

Spectra of accelerated protons and electrons 
are shown in Fig. 3. At the present epoch the 
maximum energy of protons accelerated in this 
SNR is about 140 TeV, while higher energy par- 
ticles have already left the remnant. The spec- 
tra of both electrons and protons at the reverse 
shock show significant bumps just before the cut- 
offs. Because the reverse shock presently moves in 
the rarefied medium of the ejecta, the amount of 
freshly injected low-energy particles is relatively 
small, while many high energy particles acceler- 
ated earlier have not left the reverse shock yet, 
and continue to be accelerated. 

The spectra for an early epoch, t — 100 yr are 
shown in Fig. 4. Since the shock speed was Vt — 
1.3 • 10 4 km s , the maximum energy of particles 
was also higher. For the forward shock it is about 
650 TeV. 

Spatially integrated proton and electron spec- 
tra at the present epoch t = 1620 yr are shown in 
Fig. 5. We also show the spectrum of run-away par- 
ticles. These particles have already left the simu- 
lation domain through an absorbing boundary at 
r = 2Rf. The sum of the proton spectra shown 
is the total cosmic ray spectrum produced in this 



SNR over the last 1620 years after SN explosion. 
For this SNR, the spectrum of cosmic ray protons 
have a maximum at 800 TeV. This is somewhat 
smaller than the required value to explain the knee 
in the cosmic ray spectrum. 

The results of calculations of radiation pro- 
duced by accelerated electrons and protons are 
presented and discussed in the next sections. The 
gamma-ray spectra from proton-proton interac- 
tions are ca lculated using the formalism of Kelner 
ct al. 2006. For calculations of IC gamma-rays we 
use st andar d expressions (see e.g. Blumenthal & 
Gould 1970l ). In calculations we take into account 
the target photons of the microwave background 
radiation only. We checked that other diffuse ra- 
diation fields do not contribute significantly to the 
gamma-ray production, unless the SNR is located 
in a regions with significantly (by an order of mag- 
nitude) enhanced optical and infrared radiation 
background. For calculations of synchrotron radi- 
ation we take into account that in real situations 
the magnetic field has some probability distribu- 
tion P{B). This makes the cut-off in the spectrum 
of synchrotron radiation somewhat smoother (see 
Appendix) . 

4. Hadronic scenario 

The results of calculations of the broad-band 
emission for the case of hadronic origin of high 
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Fig. 6. — The results of modeling of of nonthermal radiation of RX J1713. 7-3946 within the hadronic scenario 
of gamma-ray production. The following basic parameters are used: t = 1620 yr, D — 1.2 kpc, nn =0.09 
cm" 3 , Esn = 2.7 • 10 51 erg, M ej = 1.5M©, M f A = M b A = 23, £ = 0.05, the electron to proton ratios at the 
forward and reverse shocks K[ p = 10~ 4 and K h ep = 1.4 • 10~ 3 . The calculations lead to the following values of 
the magnetic fields and the shock speeds at the present epoch: the magnetic field downstream of the forward 
and reverse shocks Bf = 127 /jG and Bb = 21 /jG respectively, the speed of the forward shock Vf — 2760 
km s — 1 , the speed of the reverse shock Vb = —1470 km s _1 . The following radiation processes are taken 
into account: synchrotron radiation of accelerated electrons (solid curve on the left), IC emission (dashed 
line), gamma-ray emission from pion decay (solid line on the right), thermal bremsstrahlung (dotted line). 
The input of the reverse shock i s shown by the corresponding thin lines. E xperi mental data in gamma-ray 
(HESS; Aharonian et al. l2007ah and X-ray bands (Suzaku; Tanaka et al. 120081 ). as well as the radio flux 
22 ± 2 Jy at 1.4GHz (ATCA; Acero al. [2009) from the whole remnant are also shown. 



energy gamma-rays are shown in Fig. 6. The prin- 
cipal model parameters used in calculations are 
described in the figure caption. Note that at the 
present epoch already 70 % of the explosion energy 
of 2.7 x 10 51 erg has been transferred to acceler- 
ated protons (see Fig.l), most of them still con- 
fined in the shell of the remnant (see Fig. 5). The 
maximum energies are 140 TeV and 23 TeV for 
protons and electrons respectively. Although the 
reverse shock contributes significantly to the high- 
est energy particles, especially around 100 TeV 
(see Fig. 3), the gamma-ray production related to 
this population of protons is suppressed because 
of the low density of the ejecta plasma. The 
contribution of electrons directly accelerated by 
the reverse shock to synchrotron radiation and IC 
gamma-rays is more significant. In particular, the 
reverse shock produces 16% of X-rays. Moreover, 
the bump at the end of IC spectrum is contributed 



mostly by these electrons. The injection efficiency 
of electrons is adjusted in order to reproduce the 
total intensity of synchrotron X-rays. The electron 
injection efficiency at the reverse shock was ad- 
justed to reproduce the observable radio-intensity 
of the inner rin g. Th e gray scale radio i mages 



of Ellison et al. 2001 and Lazcndic et al. 2004 



were used in order to obtain a rough estimate 
of radio-flux 4 Jy at 1.4 GHz from this region. 
The modelled total radio flux is slightly below 
the observed flux. Note however that the space- 
integrated faint radio flux contains uncertainties, 
in particular the contribution of the background 
thermal radio emission can be quite significant. 

In Fig. 6 we show also the energy flux of the 
thermal bremsstrahlung. It has a maximum at 0.4 
keV and approximately equals to the energy flux 
of gamma-rays produced in proton-proton interac- 
tions . Note that the ratio of the energy flux of the 
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Fig. 7. — The same as in Fig. 6, but for the supernova shock propagating in the medium with r 2 density 
profile, assuming the following parameters t = 1620 yr, D = 1.3 kpc, E$n = 2.5 • 10 51 erg, M e j — 1.4M Q , 
M{ = M\ = 23, Co = 0.05, K[ p = 1.5 • 10~ 4 , K b ep = 4.4 • 10~ 4 . The undisturbed hydrogen number density at 
the present shock position is uh =0.12cm -3 . The calculations lead to the following values of the magnetic 
fields and the shock speeds at the present epoch: the magnetic field downstream of the forward and reverse 
shocks Bf — 131 /iG and = 28 /iG, respectively, the speed of the forward shock Vt = 2260 km s _1 , the 
speed of the reverse shock Vt, — —3010 km s _1 . 



thermal emission to the flux of " hadronic" gamma- 
rays does not depend on the density of ambient 
plasma, but depends on the shock speed and the 
cosmic ray acceleration efficiency (see Appendix). 

It is expected that a circumstellar medium 
around progenitors of Ib/c and lib supernova is 
strongly nonuniform. At the main sequence (MS) 
phase the stellar wind of progenitor creates rar- 
efied bubble in the surrounding medium. Later the 
part of this bubble is filled with a dense gas ejected 
by progenitor at the Red Super Giant (RSG) phase 
of the stellar evolution. At the end of the progeni- 
tor evolution this gas most likely will be swept up 
and transported to the MS bubble periphery by a 
powe rful Wolf-Rayet stellar wind, (see Chevalier 
120051 for a review) . 

It seems that the corresponding MS bubble 
is indeed found. SNR RX J1713. 7-3946 is sur- 
rounded by a massive shell of molecular gas (Fukui 
2008) . The densest cores of molecular gas (clouds 



C and D) are probably swept up by the forward 
shock of the SNR. This is the reason to consider 
the SNR evolution in the circumstellar medium 
with a positive density gradient. The results ob- 



tained for the gas density profile proportional to 
the square of the radius are shown in Fig. 7. A 
similar "bubble" model, but with a sharper den- 
sity p rofile was considered by Berezhko & Volk 



2006 



Results of calculations of the broad-band emis- 
sion presented in Fig. 7 are quite similar to the ones 
shown in Fig. 6, however the agreement with obser- 
vations is achieved with different set of model pa- 
rameters. Also, note that in this case the reverse 
shock produces 35% of the total X-ray emission. 

Note that in both cases the results are obtained 
for a low density gas. A similar number density 
close to i%h ~ 0.1 c m" 3 has been suggested by 
Morlino et al. 120081 in their hadronic models of 
gamma-emission of SNR RX J1713.7-3946. The 
low density gas is required also from the upper 
limit on the thermal X-ray emission from this su- 
perno va remnant (see e.g. Cassam-Chena'i et al. 



2004 Takahashi et al. I2008L Tanaka et al. 2008) 



and is in agreement with a simple estimate of the 
gamma-ray intensity (see Appendix C). Signifi- 
cantly denser gas, nn~ 1 cm" 3 , was assumed 
by Berezhko & Volk l2006i This is because they 
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Fig. 8. — Broad-band emission of RX J1713. 7-3946 for the leptonic scenario of gamma-rays with a non- 
modified forward shock. The principal model parameters are: t = 1620 yr, D = 1.5 kpc, njj =0.02 cm~ 3 , 
E SN = 1.2 • 10 51 erg, M ej = 0.74M©, M f A = 69, M\ = 10, £ = 0.1, K[ p = 2.3 • 10" 2 , K h ep = 9 • 10" 4 . The 
calculations lead to the following values of the magnetic fields and the shock speeds at the present epoch: 
the magnetic field downstream of the forward and reverse shocks Bf = 17 /iG and Bb = 31 /iG, respectively, 
the speed of the forward shock Vf — 3830 km s _1 , the speed of the reverse shock V& = —1220 km s _1 . The 
following radiation processes are taken into account: synchrotron radiation of accelerated electrons (solid 
curve on the left), IC emission (dashed line), gamma-ray emission from pion decay (solid line on the right), 
thermal bremsstrahlung (dotted line). The input of the reverse shock is shown by the corresponding thin 
lines. 



assumed that the cosmic ray acceleration occurs 
only at 1/5 part of the shock surface. 

5. Leptonic scenario 

The observed X-ray to gamma-ray energy flux 
ratio close to 15 determines the average value of 
the magnetic field B ~ 12 /iG when the leptonic 
origin of gamma-emission in RX J1713. 7-3946 be- 
comes possible. How realistic is such a week mag- 
netic field? 

The amplified magnetic fields can be quite 
weak if the SNR shock propagates in the rar- 
efied medium. However since the magnetic field 
at the reverse shock is several times lower than the 
magnetic field at the forward shock (see previous 
Section) and the X-ray intensities of the shocks 
are comparable, the IC emission will be mainly 
produced at the reverse shock; this seems to be in 
conflict with HESS observations. 

For this reason we consider a unmodified for- 



ward shock of SNR. If the injection is not effective 
at the forward shock, the energy density of cos- 
mic rays will be low and the same will be true 
for the amplified magnetic field. Such a situa- 
tion is possible, in particular, for a perpendicular 
SNR shock propagating in the medium with an 
azimuthal magnetic field. In this case we use a 
very small injection parameter i]f — 10~~ 5 at the 
forward shock and the standard value rjt, — 10~ 2 
for the injection parameter at the reverse shock. 
Thus in this scenario we deal with unmodified for- 
ward shock and modified reverse shock. Since, un- 
der these assumption, the electron to proton ratio 
K ep is close to the observed ratio for galactic cos- 
mic rays, 0.01 at 10 GeV, the electron injection 
was taken to be independent on the shock veloc- 
ity. We also use a higher value MjJ = 69 for the 
forward shock, and adopt the distance D — 1.5 
kpc that is close to an up per limit a ccord ing to 
Cassam-Chenai' et al. 120041 and Fukui l2008l 

We should note that the shock speed is one of 
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Fig. 9. — The same as in Fig. 8 but for unmodified forward and reverse shocks. The principal model 
parameters are: t = 1620 yr, D = 1.5 kpc, n H =0.02 cur 3 , E SN = 1-25 • 10 51 erg, M ej = 0.82M©, M b A = 23, 
M{ = 69, Co = 0.1, K[ p = 2.0 • 10~ 2 , K b ep = 6.8 • 10~ 2 . The calculations lead to the following values of the 
magnetic fields and the shock speeds at the present epoch: the magnetic field downstream of the forward 
and reverse shocks Bf = 17.5 and Bf, = 13 fiG, respectively, the speed of the forward shock Vf = 3870 
km s _1 , the speed of the reverse shock Vf, = —1290 km s _1 . 



the key parameters in order to achieve the good 
feet. The shock speed is related with a remnant 
size at given age of the remnant. This explains a 
slight difference of the assumed value of distances. 

The results are shown in Fig. 8. The magnetic 
fields at the reverse shock and the forward shock 
are comparable in this case. The energy of cosmic 
rays is only 10% of the explosion energy. The cos- 
mic ray pressure at the forward shock is only 2% of 
the ram pressure poR 2 f- The particles are mainly 
accelerated at the reverse shock. The maximum 
energies are 55 TeV and 36 TeV for protons and 
electrons respectively. The reverse shock produces 
16 % of X-ray emission. 

For completeness we also consider the case 
when both forward and reverse shocks are not 
modified. The injection efficiency T]b = r)f = 10~ 5 
is used. The results are shown in Fig. 9. The 
energy of cosmic rays is 3% of the explosion en- 
ergy. The reverse shock produces 20% of X-ray 
emission. 
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Fig. 10. — Radial profiles of 2 keV X-rays (multi- 
plied to 0.04, solid line), 1 TeV gamma-ray emis- 
sion (dashed line) and 1.4 GHz radio-emission 
(multiplied to 10 3 , dotted line), calculated for the 
hadronic scenario in the uniform medium. 
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6. Radial profiles 



The radial profiles of brightness distributions 
of X-ray, gamma-ray and radio- emissions calcu- 
lated for leptonic and hadronic models are shown 
in Fig. 10-12. The projection effect is taken into 
account. 

For the hadronic scenario in the uniform 
medium, all three components of radiation in the 
radio, X-ray and gamma-ray bands peak at the 
forward shock (approximately 0.5° of the angular 
radius. Note that the synchrotron emission has 
a sharper feature at the forward shock than the 
gamma-ray emission. The synchrotron emission 
also shows a noticeable peak at the reverse shock 
(approximately 1/3 of the brightness at the for- 
ward shock). The slight feature of gamma-rays 
(at the level of 10 %) at the reverse shock is due 
to the IC gamma-rays. The middle peak of X-ray 
emission in Fig. 11 is explained by the electrons 
accelerated at the forward shock and reaching the 
contact discontinuity with a stronger magnetic 
field. In Fig. 13 the calculated gamma-ray pro- 
files are compared with the HESS data. 

The most prominent feature at all X -ray radial 
profiles from Figs. 10-12 is a clear visible filament 
at the reverse shock. Note that the projection ef- 
fect works differently for reverse shock in compar- 
ison with the forward shock (see Appendix D). As 
a result, at the reverse shock, the filaments are 
visible even in the relatively weak magnetic fields 
(13-30 fj,G). Such filaments have been indeed ob- 
ser ved b y Chandra in th e inne r rim (Uchiyama et 
al. I2003L Lazendic et al. l2004h . 

The scale Id of an exponential decrease of the 
synchrotron emissivity downstream of the shock 
with a veloc ity V a is given by (Zirakashvili & Aha- 
ronian[2007|) 
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Here k = B u /Bd is the ratio of magnetic fields 
upstream and downstream of the shock. The only 
available profile of the linear fila ment i n the inner 
region is given by Lazendic et al. |2004 Note that 
the profile is probably slightly influenced by CCD 
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Fig. 11. — Profiles of 2 keV X-ray emission (mul- 
tiplied to 0.1, solid line), 1 TeV gamma-emission 
(dashed line) and 1.4 GHz radio-emission (multi- 
plied to 10 3 , dotted line) for the leptonic scenario 
with the non-modified forward shock. 
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Fig. 12. — Radial profiles of 2 keV X-rays (multi- 
plied to 0.1, solid line), 1 TeV gamma-ray emission 
(dashed line) and 1.4 GHz radio-emission (multi- 
plied to 10 3 , dotted line) for the leptonic scenario 
with the non- modified forward and reverse shocks. 
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Fig. 13. — Radial Profiles of 1 TeV gamma-rays 
for the hadronic scenario in the uniform medium 
(solid) and for the leptonic scenario with the un- 
modified forward shock (dashed). The profiles 
smoothed with a Gaussian point spread func- 
tion with a ~ 0.05° are also shown (thin lines). 
The triangles show the azimuthally averaged TeV 
gamma-ray r adial p rofile observed by HESS (Aha- 
ronian et al. [2007a). 



gap of Chandra. The filament width at the half 
of magnitude is about 40". This corresponds to 
30" for the exponential scale (see Appendix D) or 
Id = 7 • 10 17 cm at D = 1.5 kpc. Using the last 
formula we find Bd = 25 fiG . So the magnetic 
field at the reverse shock in our calculations corre- 
sponds to the width of X-ray filaments. According 
to Fig. 10, filaments with a similar width must be 
also observed at the outer X-rims in the hadronic 
model. 



7. 



Maximum energies and magnetic am- 
plification at the reverse and forward 
shocks of SNR 



According to Zirakashvili k Ptuskin (|2008f) the 
maximum energy of protons E m accelerated at the 
quasi-parallel shock moving in the circumstellar 
medium with the magnetic field strength Bq and 
Alfven velocity V a — Bq/^/Attpq may be found 
from the expression 



Err, = 21 TeV 
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Here r\ esc = 2Fe/poVj is the ratio of the energy 
flux of run-away particles Fe to the flux of kinetic 
energy p Vj /2, Bb is the strength of the random 
magnetic field in the circumstellar medium, u* = 
1.73 • 10 3 km s _1 . The expansion parameter m eX p 
is determined as m cxp = din R t / d In i and is close 
to 0.5 in supernova remnants in the transition to 
the Sedov phase in the homogeneous medium. The 
speed u' is a so-called normalized shock velocity 
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Eq. (12) is derived for the diffusive shock accel- 
eration in the presence of the non-res onant stream- 
ing instability suggested by Bell (|2004r ). This 
equation is valid for u' > u*. The maximum en- 
ergy for lower velocities when the resonant insta- 
bility should be taken into account is given by the 
same expression with the logarithm in the denom- 
inator only. The ratio rj esc = 0.14 for the cosmic 
ray modified shock with the compression ratio a = 
6 and the compression ratio of the thermal sub- 
shock 2.5. For magnetic field strength Bq = 5/iG 
in the circumstellar medium, Vf = 2760 km s _1 
from the hadronic model (see Fig. 6) and for a char- 
acteristic value of In 2Bq /B = 5 we find u' = 1930 
km s _1 and E m = 136 TeV. This normalized shock 
velocity also determines the strength of the ampli- 
fied magnetic field B ~ 2.5B = 12. 5/iG upstream 
of the forward shock or B ~ 75/iG downstream 
of the shock after the compression in the shock 
transition region. 

These numbers are in a qualitative agree- 
ment with our numerical results obtained for the 
hadronic scenario. We found using Eq. (12) that 
the maximum energies are significantly (a factor 
of 40) lower in the leptonic scenario with the non- 
modified forward shock. This is because a small 
amount of accelerated protons. However the max- 
imum energy is underestimated if one uses Eq. 
(12) for an oblique forward shock. The particles 
that run-away from such a shock will produce a 
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stronger electric current in a comparison with the 
quasi-parallel case. Since the maximum energy 
E m is proportional to the current which drives 
the streaming instability this will permit to reach 
the maximum energy of several tens TeV for the 
leptonic scenario. This makes the leptonic sce- 
nario with the non-modified forward shock to be 
self-consistent. 

The maximum energy and the strength of the 
amplified magnetic field at the reverse shock de- 
pend on the magnetic field of ejecta. The regular 
magnetic field strength B ~ 100 G of the com- 
pact progenitor of Ib/c supernova with a radius 
R = 10 12 cm will drop down to 10~ 12 G after the 
homogeneous radial expansion of the ejecta up to 
radius 3 pc~ 10 19 cm. Probably inside the progen- 
itor there exist stronger random magnetic fields of 
the order B ~ 10 4 G. This fields will drop down to 
10 -10 G after the homogeneous expansion. In the 
real situation the expansion is not homogeneous 
and the random magnetic field is stretched in the 
radial direction. This may additionally amplify 
the field by a factor of ~ 10 that is the ratio of 
the stellar radius to the expected azimuthal inho- 
mogeneity scale of the velocity field. This gives an 
upper limit 10~ 9 G for the magnetic field strength 
of ejecta. The radial orientation of the field is fa- 
vorable for the ion injection at the reverse shock. 

The close value one may find in the model 
of magneto-rotational super nova explosion of 
Bisnovatyi-Kogan et al. 2008. The magnetic field 
is amplified in the differentially rotating iron core 
up to B ~ 10 17 G after the core collapse and 
pushes the core material outward producing the 
supernova explosion. For the core radius 10 6 cm 
we found B ~ 10~ 9 G after the homogeneous ex- 
pansion and B ~ 10 -8 G after inhomogeneous 
expansion. We shall use this upper optimistic 
limit Bq ~ 10~ 8 G for the estimate below. 

In the hadronic model the reverse shock moves 
with the speed u — —4500 km s" 1 in the ejecta 
frame of reference. The ejecta number density 
n< H = ^ ' 10~ 4 cm~ 3 . Then the Alfven veloc- 
ity V a — 0.67 km s^ 1 and the normalized shock 
velocity v! — 5.7 • 10 9 km s _1 for rj esc = 0.14. 
Then we find E m = 2.3 TeV using Eq. (12). The 
magnetic amplification fa ctor w as not determined 
by Zirakashvili & Ptuskin[2008 for such high nor- 
malized velocities. The amplification factor is 190 
for the normalized velocity u' = 3.9 • 10 9 km s _1 . 



However we find B ~ 300£?o = 3/iG after an ex- 
trapolation. So the amplified field downstream of 
the shock is B ~ 18/iG. The maximum energy 
and amplified magnetic field may be even higher 
if one takes into account that the number density 
of the highest energy particles is higher in com- 
parison with what is expected at the shock with 
run-away particles. These particles cannot run- 
away from the reverse shock easily because the 
upstream region of the reverse shock is located in- 
side the remnant (see also bumps at the end of the 
spectra in Fig. 3). Amplified magnetic field weakly 
depends on the initial seed field Bq as B cx B^ 4 
while the dependence of the maximum energy is 
stronger E m oc B ^j 4 for small Bq (see Zirakashvili 



& Ptuskin[200ifor details). 

We conclude that for the magnetic field strength 
Bq ~ 10 _8 G of ejecta the reverse shock may accel- 
erate particles up to multi-TeV energies and may 
amplify magnetic fields up to tens of fiG in the 
remnant considered. 

8. Discussion 

8.1. Thermal X-ray emission 

The nonthermal emission of SNRs cannot be 
treated out of the context of thermal emission. 
In particular, any standard hadronic model of 
gamma-rays predicts heating of the gas to high 
temperatures with intense X-ray emission. The 
lack of such emission from the most pronounced 
TeV emitting SNR, RX J1713. 7-3946, can be inter- 
preted as an argument against the ha droni c mod- 
els of gamma-rays (Katz & Waxman 120081 ) or an 
evidence of extremely high efficiency of transfor- 
mation of the kinetic energy of explosion to accel- 
eration of relativistic particles (Drury et al. 2008) . 

Though our numerical calculation show that in 
the hadronic models of gamma-rays the thermal 
bremsstrahlung appears to be an order of magni- 
tude below the Suzaku X-ray data (see Figs 6 and 
7) , the actual flux of X-ray emission in this energy 
range is higher due to the contribution fr om th e 
line X-ray emission (see e.g. Ellison et al. 120071 ). 
The latter strongly depends on the chemical com- 
position of the X-ray emitting plasma. For the 
composition similar to one observed in the Sun, 
the intensity of X-ray lines can significantly excee d 
the free-free continuum (see Ellison et al. 120071) . 
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This means that such lines must be observed con- 
trary to X-ra y data . Based on this fact, Cassam- 
Chenai et al. (20q3) derived an upper limit on the 
ambient plasma density uh ~ 0.02 cm -3 . This 
value is several times lower than nu ~ 0.1 cm -3 
in the hadronic models considered in this paper. A 
simplest explanation of this could be that chemi- 
cal composition of the circumstellar gas is strongly 
different from the solar one e.g. due to absorption 
of X-ray line emitting ions by dust grains. 

We should note that the thermal X-ray emis- 
sion may be strongly depressed at the supernova 
shocks that are very strongly modified by the cos- 
mic ray pressure when the shock compression ratio 
a >> 10. In this case only a small part of dissi- 
pating energy is transferred to the thermal plasma 
and the gas temperature may be below 0.1 keV 
even under t he ele ctron-ion thermal equilibrium 
(Drury et al. 2008) However such compression ra- 
tios are not observed in our model. The strong 
heating of the upstream plasma due to the damp- 
ing of AlfVen waves results in the modest shock 
modification with a ~ 6 — 7. Even if one assumes 
that the waves are not damped and the gas heat- 
ing is negligible, the wave pressure will be high 
and it will also prevent t he str ong shock modifi- 
cation (see Caprioli et al. 12003 for details). Very 
high compression ratios are possible at radiative 
shocks in a high density plasma where the waves 
are damped and the gas loses energy radiatively. 

Another possibility to explain a low thermal X- 
ray luminosity of the remnant is related with a 
nonuniform circumstellar medium. Chandra ob- 
servations reveal a complex network of non-steady 
bright spots and filaments i nside the X-ray emit- 
ting shell (Uchiyama et al. 120070 . If dense fila- 
ments with a small volume factor exist already in 
the circumstellar medium and determine the mean 
density of the plasma, say 1 cm -3 , they will not be 
shocked by the SNR shock coming and the mean 
density of the shocked material will be rather low 
say 0.01 cm -3 . However the plasma density "seen" 
by TeV accelerated protons can be close to the to- 
tal mean value because of the relatively fast diffu- 
sion. As a result they will be able to produce TeV 
gamma-ray with a flux close to the observed one. 
An enhanced magnetic field inside the filaments 
may also produce the observed X-ray variability 
of bright spots (Uchiyama et al. 120071 ). 

An interesting consequence of this scenario is 



that the gamma-ray intensity at low energies can 
be reduced because of the slow diffusion of GeV 
protons which prevents their effective penetration 
into the dense filaments. As a result, the gamma- 
ray spectrum due to p-p interactions can be quite 
similar to the IC gamma-ray spectra expected in 
leptonic models. This possibility makes quite dif- 
ficult the identification of the origin of gamma- 
radiation. 

It is interesting to note that the thermal X-rays 
produced in the gas shocked by the forward shock 
are rare observed in young supernova remnants. 
For example, the thermal X-ray emission of the su- 
pernova remnants Tycho and SN1006 is produced 
by shocked ejecta that is en riched by h eavy ele- 
ments (Cassam-Chena'i et al. 120071 120081 ). In this 
regard the SNR RX J1713. 7-3946 perhaps is not a 
very special example. It is at a later evolutionary 
phase in comparison with SNRs mentioned above, 
so the ejecta is situated deeply inside and has a 
very low density (see Fig. 2). This could be the 
reason why the thermal X-ray emission of ejecta 
is not observed. On the other hand, the thermal 
emission from the downstream region of the for- 
ward shock is not detected similar to many other 
young SNRs. 

The young remnant RCW 86 is well known 
source with thermal X-ray line emission observed 
from the forward shock (see Yamaguchi et al. 
2008) . It is interesting to note that in this remnant 
the chemical abundance of heavy ions downstream 
of the forward shock seems indeed below the solar 



8.2. Energy spectra and morphology 

The comparison of Figs. 6&7 with Figs, 
shows that the spectral shape of gamma-emission 
is better reproduced in hadronic models. However 
the difference perhaps is not sufficient for a def- 
inite conclusion concerning the origin of gamma- 
radiation radiation based on the current data. The 
spectra predicted by two models deviate signifi- 
cantly at very high (> 10 TeV) and low < 10 GeV 
energies. The first energy domain can be explored 
with the next generation ground-based detectors 
with significantly improved sensitivity compared 
to HESS. In particular, the confirmation of the 
the last spectral points above 30 TeV reported by 
HESS with limited statistical significance, would 
be a strong argument in favor of the hadronic 
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origin of gamma-rays. Also, the hadronic model 
predicts gamma-ray flux below 100 GeV at the 
level of > 10~ n erg/cm 2 s, while the IC gamma- 
ray flux predicted by the leptonic model around 
1 GeV is closer to 10~ 12 erg/cm 2 s. It is ex- 
pected that the measurements conducted with the 
AGILE and Fermi telescopes will soon provide 
a decisive answer for the MeV/GeV flux from 
RXJ 1713.7-3946. The detection of a signal close 
to 10 -11 erg/cm 2 s would imply hadronic origin of 
radiation. On the other hand, the much lower flux 
would favor the IC origin of radiation but still 
could not reject the hadronic origin, given that 
gamma-rays could be produced in dense filaments. 
As discussed above, because of slow diffusion of 
low energy protons that prevents their effective 
penetration into the filaments, the GeV emission 
from p-p interactions in the dense filaments could 
be suppressed. This model can also give a rea- 
sonable, in our view, explanation of the lack of 
thermal X-ray line emission. 

We should note, however, that the real pic- 
ture could be more complex, in particular we may 
expect non-negligible deviations of energy spec- 
tra because of non-spherical geometry, different 
shock properties at different parts of the shock sur- 
face, etc. Moreover, one may expect comparable 
contributions of gamma-rays from interactions of 
both accelerated protons and electrons, but pro- 
duced in different proportions in different regions 
of the remnant (see below). In this regard, de- 
tailed spectroscopic measurements of gamma-rays 
from different parts of the remnant seem to be 
a key towards understanding of acceleration pro- 
cesses which take place in this source. 

If the magnetic field is indeed as high as re- 
quired in hadronic models, we should expect, ac- 
cording to Fig. 10, thin filaments with a width of 
about 30" in the remnant periphery. The lin- 
ear filaments are ind eed pr esent at Chand ra im- 



age (Uchiyama et al. 120031 Lazendic et al. 2004) 
but in the vicinity of the inner ring that might 
be related to the reverse shock. This agrees with 
Fig. 10-12. The most r ecent X -ray image of XMM- 
Newton (Acero et al. 120091 ) shows several linear 
filaments with the width 1 — 2' present in this im- 
age, but they all are situated inside the boundary 
of the X-ray emission. The boundary itself is not 
sharp that is presumably related with a nonuni- 
form circumstellar medium. The lack of sharp X- 



ray boundary is better explained by leptonic mod- 
els. According to Fig. 11 and Fig. 12 some amount 
of X-rays is expected from the region located be- 
yond the forward shock. In the leptonic model 
it is possible because accelerated electrons have 
sufficiently high energies to produce synchrotron 
X-rays in the low magnetic field. 

The gamma-ray radial distributions predicted 
by hadronic models are sharper than the distri- 
butions calculated within the leptonic models (see 
Fig. 13). However, because of limited angular reso- 
lution of gamma-ray telescopes, it is not possible, 
unfortunately, to distinguish between the distribu- 
tions predicted by two models. It is demonstrated 
in Fig. 13 in which the profiles are smoothed with 
a typical for Cherenkov telescope arrays point 
spread function assuming Gaussian distribution 
with a — 0.05°. Both smoothed profiles reason- 
ably agree with the angular distribution of TeV 
gamma-rays as measured by HESS. An improve- 
ment of the angular resolution of future detectors 
by a factor of two should allow identification of 
the origin of the parent particles. On the other 
hand, this can be done with current detectors like 
HESS, VERITAS and MAGIC, for a young nearby 
SNR with a larger angular size. In this regard, the 
young supernova remnant Vela Jr (RX J0852.0- 
4622) with gamma-ray properties quite similar to 
RX J 1713/7-3946 as reported by HESS (Aharonian 
et al. l2007bh . but with angular diameter w 1°, is 
a promising object for such studies. 

8.3. Shocked clouds in the shell? 

The observed ratio of the forward shock radius 
to the reverse shock radius robustly constraints 
the type of supernova explosion that produced 
SNR RX J1713. 7-3946. For reasonable values of 
the energy of supernova explosion E$n < 3 • 10 51 
erg, the ejecta mass must be small, M e j < 2M , 
otherwise the forward shock speed would be too 
low for production of synchrotron radiation ex- 
tended to hard X-ray domain. Such small ejecta 
masses correspond to lb/ c or lib core collapse su- 
pernova (Chevalier 12005 ). Circumstellar medium 
around these types of supernova is created by the 
stellar wind of the supernova progenitor. It may 
be a low density (n# ~ 0.01 cm~ 3 ) bubble cre- 
ated by stellar wind of progenitor during the main 
sequence period of the star evolution or by a stel- 
lar wind of a Wolf-Rayet progenitor. Because of 
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Fig. 14. — Broad-band emission of RX J1713. 7-3946 for the composite scenario of gamma-rays with a non- 
modified forward shock and dense clouds. The principal model parameters are: t = 1620 yr, D = 1.5 kpc, 
n H =0.02 cm" 3 , E SN = 1.2 • 10 51 erg, M ej = 0.74M©, M f A = 55, M\ = 10, Co = 0.1, K{ p = 1.4 • 10~ 2 , 
K\ v = 9 • 10~ 4 . The calculations lead to the following values of the magnetic fields and the shock speeds 
at the present epoch: the magnetic field downstream of the forward and reverse shocks Bf — 22 /iG and 
Bi, = 31 [J.G, respectively, the speed of the forward shock Vf = 3830 km s _1 , the speed of the reverse shock 
Vb = —1220 km s _1 . The following radiation processes are taken into account: synchrotron radiation of 
accelerated electrons (solid curve on the left), thermal bremsstrahlung (dotted line), IC gamma-ray emission 
of the entire remnant including forward and reverse shocks (dashed line), hadronic component of gamma-rays 
from the remnant's shell (solid line on the right) as well as from dense clouds assuming the factor of 120 
enhancement of the flux (thin dashed line). We also show the total gamma-ray emission from the entire 
remnant including the dense clouds (thin solid line). 



very low densities, the observed gamma-rays can 
be produced only via IC scattering of electrons. If 
the bubble is created during the Red Supergiant 
stage of the supernova progenitor or during the 
interaction of the slow Red Supergiant wind with 
the fast Wolf-Rayet wind, the gas density could be 
sufficiently high for effective production of gamma- 
rays at interactions of accelerated protons. 

In all models considered the forward shock of 
SNR have swept only 6-15M Q . This mass is com- 
parable or slightly less than the mass ejected by 
progenitors of Ib/c and lib supernova during the 
stellar evolution. This means that the shock just 
have swept the progenitor's material ejected dur- 
ing the stellar evolution. The interaction with 
the molecular gas surrounding the remnant only 
starts. The exceptions could be the very dense 
cores of molecular gas (clouds C and D) that are 
situated already inside the forward shock (Fukui 



2008). Probably the forward shock has enveloped 
the clouds, and a high pressure gas from the down- 
stream region drove secon dary shocks into the 
clouds (see Chevalier 19771 for a review). These 
shocks are strongly decelerated and become radia- 
tive. That is why they do not produce th e ther - 
mal X-rays (see also Cassam-Chenai et al. 120041) . 



The nonthermal X-rays can be produced at these 
shocks by high-energy electrons from the rem- 
nant's shell in the compressed magnetic field of the 
cloud, while the highest energy gamma-rays from 
pion decay are mainly produced inside the c loud 
where the target density is high (Fukui [2008(). 

If the highest energy protons freely pene- 
trate into the clouds, the corresponding gamma- 
emission from the pion decay may exceed the 
gamma-emission from the pion decay of the rem- 
nant by a factor that is the ratio of the cloud mass 
to the mass swept up by the forward shock. If the 
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clouds C and D in deed contain about 100-1000 so- 
lar masses (Fukuil2008|), this ratio should be close 
to 10-100. Since this is not the case, as it follows 
from the HESS results, the hadronic models can 
survive only if the penetration into the clouds is 
depressed. Such a possibility seems unlikely, at 
least for the highest energy protons. 

On the other hand, in the leptonic models of 
gamma-rays, in which the bulk of the gamma- 
ray emission is produced via IC scattering of elec- 
trons, one may expect significant contribution of 
hadronic gamma-rays produced in the clouds C 
and D. Namely the gamma-ray fluxes from p-p 
interactions shown in Figs.8&9 can be enhanced 
by 1 or 2 orders of magnitude. Apparently, from 
these regions we should expect also low energy, 
MeV/GeV gamma-rays. Because of low density 
of the ambient gas, the gamma-radiation outside 
the C and D clouds will be dominated by the 
IC component, therefore the gamma-ray flux at 
MeV/GeV energies is expected to be faint (see 
Figs.8.&9). 

In Fig. 14 we demonstrate the feasibility of the 
"composite" scenario. The hadronic component of 
gamma-rays from dense clouds is shown by multi- 
plying the flux expected from the low density gas 
downstream of the forward shock by a factor of 
120. The ratio of the gamma-rays from clouds to 
the flux from other parts of the shell depends on 
the ratio of the mass of clouds to the mass of the 
shell, provided that all particles freely enter the 
dense clouds. This however could not be the case 
especially for the low energy particles which intro- 
duces large uncertainties in the predictions of the 
flux of this component of radiation. On the other 
hand, if the enhancement of hadronic gamma-rays 
in clouds is indeed very large, this could help to 
improve the fit of the gamma-ray data both at 
low (sub TeV) and high energy (multi-TeV) re- 
gions (see Fig. 14). Note that in this scenario we 
expect also gamma-rays from dense clouds due to 
bremsstrahlung of electrons. However, as long as 
the total energy in protons significantly exceeds 
the energy in relativistic electrons, the interaction 
of latters with the dense gas will dominate over 
the bremsstrahlung. 



8.4. Reverse shock in hadronic and lep- 
tonic models 

The parameters characterizing acceleration of 
particles in the reverse shock in the leptonic mod- 
els are not strongly constrained. While in the 
hadronic models the reverse shock must be an ef- 
ficient accelerator of multi TeV electrons in order 
to produce the synchrotron X-radiation observed 
from the inner ring, in the leptonic models the 
acceleration by the inverse shock is not required. 
Indeed, in the leptonic models the electrons ac- 
celerated at the forward shock may easily reach 
the contact discontinuity downstream of the for- 
ward shock with relatively weak magnetic fields. 
They will produce X-rays in the amplified mag- 
netic field between the contact discontinuity and 
the reverse shock (see Fig. 11). The magnetic field 
cannot be weak in this region since the inner ring is 
clearly visible in radio. The corresponding radio- 
electrons must be accelerated at the reverse shock. 
The magnetic field amplified at the reverse shock 
has probably rather small scales. This is because 
the scale of the most unstable MHD mode is pro- 
portional to the strength of the undisturbed field 
in the ca se of t he non-resonant streaming instabil- 
ity (Bell [2002). 

The magnetic field of ejecta B < 10~ 8 G is 
significantly smaller than the field in the circum- 
stcllar medium (presumably about 10 -6 G). The 
small-scale magnetic field is not a very efficient 
scatter of highest energy electrons. As a result, 
the diffusion may be significantly faster than the 
Bohm diffusion, and the acceleration at the re- 
verse shock up to TeV energies might be impossi- 
ble. However, acceleration of GeV electrons pro- 
ducing synchrotron radio emission is quite prob- 
able. On the other hand, if the reverse shock is 
an effective accelerator, it may produce the addi- 
tional component(s) of nonthermal radiation. For 
example, in the case of weak magnetic fields at 
the reverse shock, the hadronic models may re- 
sult in detectable IC gamma-rays from the central 
region of the remnant. Probably the situation is 
different in the remnant Cas A, where the mag- 
netic field at the reverse shock is stronger than 
the magnetic field at the forward shock. This may 
be related with a different acceleration efficiency 
of the shocks in the remnant Cas A. For example, 
the reverse shock may be modified and the forward 
shock not. The reverse shock energetics in Cas A 
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Fig. 15. — The synchrotron emissivity of a single 
electron in the isotropic magnetic field. The func- 
tion R(x) for single value of the magnetic field 
strength given by Eq.(A2) (dashed line) and the 
function R\(x) (solid line) given by Eq. (A5) are 
shown. 

is of the order of 10 % from the total one. It is 
enough to explain the observable TeV gamma-rays 
in the framework of hadronic models. 

9. Conclusion 

Our study of the nonlinear diffusive shock ac- 
celeration in RX J1713. 7-3946, and calculations of 
the related electromagnetic radiation lead to the 
following conclusions: 

1. The present reverse shock position may be 
reproduced if the SNR RX J1713. 7-3946 is the 
remnant of Ib/c or lib type supernova with the 
ejected mass M e j < 2M Q . 

2. The available radio, X-ray and gamma-ray 
data can be explained by both the leptonic and 
hadronic models, although with certain caveats. 
The spectral shape of gamma-emission is better 
reproduced in hadronic models. The morphol- 
ogy of electromagnetic emission in X-ray and TeV 
gamma-ray bands can be in principle reproduced 
by both the hadronic and leptonic models. 

3. The hadronic origin is possible only if the 
following conditions are satisfied: 

a) The forward shock is strongly modified through- 
out all surface. 

b) Chemical abundance of heavy ions in X-ray 
emitting regions should be significantly below the 




1 10 100 

Fig. 16. — Spectra of synchrotron X-rays pro- 
duced at the shock with a distributed magnetic 
field (Eq. (B5), solid line) and at the shock with 
a single- valued magnetic field (Eq.(B3), dashed 
line). Suzaku data (circles) for the remnant RX 
J1713. 7-3946 are also shown. 

solar one. 

c) The lack of the predicted by hadronic models 
sharp outer X-ray rim can be explained by the 
nonuniform circumstellar medium. 

d) If the forward shock interacts with clouds D and 
C, the penetration of the highest energy protons 
into the clouds must be suppressed. 

4. Any detection of GeV gamma-rays by AG- 
ILE and Fermi satellites at the level of F > 10~ n 
erg cm~ 2 s , especially outside the dense clouds 
C and D, would be a decisive argument in favor 
of the hadronic origin of TeV gamma-rays. The 
leptonic models predict very hard energy spectra 
of gamma-rays with absolute energy fluxes as low 
as F = 2 • 10~ 12 erg cm~ 2 s -1 . However a low flux 
and hard gamma-ray spectrum are expected also 
in the hadronic models with dense filaments, as- 
suming energy-dependent cosmic ray penetration 
into the clouds. 

5. The leptonic models imply that the bulk 
of the observed gamma-radiation is produced via 
IC scattering of electrons. However, these mod- 
els allow also significant contributions of hadronic 
gamma-rays from dense gas region in the shell, in 
particular from the clouds C and D. In such "com- 
posite" models we should expect relatively high 
fluxes of low energy gamma-rays from the C and 
D clouds, while outside these condensations the 
MeV/GeV flux is expected to be quite low. 
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Fig. 17. — The function b(x) that describes the ra- 
dial profile of the X-ray brightness for the forward 
shock (Eq. (D2), solid line) and for the reverse 
shock (Eq.(D3), dashed line). 



6. If the inner ring of RX J1713. 7-3946 is indeed 
produced by the reverse shock, its proper motion 
can be measured. X-ray data of Chandra may be 
used for this purpose since the X-ray filaments in 
the inner ring region have bee n obs erved (Lazendic 
et al. I2004L Uchiyama et al. l2003h . In all models 
considered the reverse shock moves with several 
thousand kilometers per second to the center of 
the remnant. The model with r 2 density profile 
has the highest absolute value of the reverse shock 
speed, Vb — —3010 km s -1 or even larger for a 
sharper density profile. The measurement of the 
reverse shock speed will constrain the gas density 
distribution around this SNR. 
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A. Synchrotron emission in the turbulent medium 



The synchrotron emissivity in the medium with nonuniform magnetic field can be written as 

e(u) = / dBBP(B) [ p 2 dpN{p)R{oj/u c ) (Al) 



where to c is the characteristic frequency of synchrotron radiation lo c — 1.5 qBp 2 /m^c 3 . The function P(B) is 
the probability distribution of the magnetic field. If the magnetic field strength has a single value Bo, then 
P{B) = 5{B — Bq) and this equation is reduced to a standard formula for the synchrotron emissivity. 

The function R(x) of the argument x = uj/lu c describes synch rotron radiation of a single electron in 



magnetic field with chaotic directions. Crusius & Schlickeiser (1986) derived an exact expression for R(lu/lo c ) 
in terms of Whittaker's function. With an accur acy of several percent R(uj/uj c ) can be presented in a simple 
analytical form (Zirakashvili & Aharonian 120071 ) : 



= l-Slcxpt-s) 
Vx- 2 /3 + (3.62A) 2 

In the turbulent medium the magnetic field is not uniform and is described by a probability distribu- 
tion P(B). It reve als th e exponential tails when the magnetic field is amplified by turbulent motions (see 
Schekochihin ct al. 2004). The probability distribution found in the numerical simulations o f the n onresonant 



streaming instability may be written in the simple analytical form (Zirakashvili & Ptuskin 120081 ): 

P(B) = -J^- exp(~V6B/B rms ). (A3) 

rms 

1 /2 

Here B rms = (£> 2 ) is the square root of the mean square of the random magnetic field (note the 
corrected missprint in the normalization factor in this equation) . Now the synchrotron emission of a single 
electron is described by the function R\{x): 

oo 

Rx{x) = 6 J dyy 2 R(x/y)eM-^y)- (A4) 
o 

We found a useful analytical approximation of this function 

R^x) = 1.50a; 1 / 3 (l + 1.53x 1/2 ) exp (-^^x 1 ' 2 ^ (A5) 
The functions R(x) and Ri(x) are shown in Fig. (15). 

Qualitatively similar results were found by Bykov et al. 120081 but assuming Gaussian probability distri- 
bution for P(B). 



B. Analitical expressions for the synchrotron X-rays 

The analytical approximations f or th e electron and synchrotron spectra at the non-modified shock were 
found by Zirakashvili & Aharonian l2007l . The synchrotron loss dominated regime was considered. In the case 
of the single valued magnetic field the synchrotron spectra are parameterized with the principal parameter 

81 fmiecVf _ 2.2 keV / Vf 
~ 64(1 + kV2)2 q 2 VB - (1 + ^1/2)2^ (^.103 kms -i 
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Here k is the ratio of the upstream and downstream magnetic field strength. In the case when magnetic field 
has the same value upstream and downstream of the shock (k = 1) the spectrum of synchrotron X-rays has 
the following form 



F(E) cx E~ 



1 + 0.46 



0.6 



11/4.8 



exp 




(B2) 



In the case when ratio of the upstream and downstream magnetic field strength k = 1/vTT, the spectrum 
of synchrotron X-rays is mainly determined by the downstream region and has the following form 



F(E) cx E~ 



0.5 



11/4 




CX P ~\ TT- 



(B3) 



Suzaku data obtained for RX J1713 . 7-3946 in a br oad en ergy interval from 0.4 to 40 keV are well fitted by 
these spectra (Uchiyama et al. 120071 Tanaka et al. l2008h . For example, E sa 0.94 keV for Eq.(B3). 

In the case of distributed magnet ic fiel d one should use Eq.(A5) instead of Eq.(A2). Using the electron 
spectra of Zirakashvili & Aharonian 120071 we found the following synchrotron spectra 



F{E) cx E~ 



E \ 46 ' 
l+ n.I72(-j 



25/(12-0.46) 



exp 



and 



F(E) cx E~ 2 
Here the parameter E\ is 

F x = 



1 + 0.185 



V6 



0.4-1 25/(12-0.4) 



hm e cVj 



exp 



0.134 keV 



E 
~E~i 

1/3 



1/3' 



K = 1, 



, K = 1/V11 



128(1 + n i / 2 y 



(1 + kV2)2^ b V3-10 3 kms 



V f 



(B4) 



(B5) 



(B6) 



Suzaku data may be fitted with E\ as 0.036 keV using the formula (B5). This value and t\b = 2 corresponds 
to the shock velocity Vf — 3300 km s _1 . This number is slightly above the corresponding values in the 
hadronic scenario when the larger compression ratio of the modified shocks and the additional input of the 
faster reverse shock should be taken into account. The shock velocity is higher in the leptonic scenario when 
the adiabatic losses downstream of the shock play a role. 

The spectra given by Eqs.(B3),(B5) and Suzaku data are shown in Fig. 16. The distributed magnetic 
field makes the cut-off of synchrotron spectra slightly smoother in comparison with a single-value case. The 
cut-off energy is a factor of 1.6 higher for the distributed magnetic field. 

The spectrum (B5) shows a rather slow drop beyond the cut-off. It may be used for explanation of the 
high-energy X-ray tail observed in the supernova remnant Cas A. 



C. Simple estimates of the gamma-ray emission and thermal X-ray emission in SNR. 

Using Eq. (10) one can estimate the characteristic electron temperature downstream of the forward shock. 
Assuming a linear radial profile of the gas velocity we found the following expression for the maximum electron 
temperature 
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(CI) 

Here £ ff = P ff / Po^/ is the ratio of the gas pressure downstream the shock P g to the dynamical pressure poV 2 
and a is the total compression ratio of the shock. The time derivative of the temperature was neglected in 
Eq.(lO). For cosmic ray modified shock in the hadronic model £ g ~ 0.25, a — 6.5 and we obtain T e = 0.6 
keV in agreement with the numerical results. 

It is clear from Eq.(Cl) that the dependence of the electron temperature on the shock parameters is 
rather weak. E.g at the Sedov stage T e oc t~°' m . It also weakly (only a factor of 2) affected by the shock 
modification. Note that the electron-ion equilibration case when the electron temperature T e oc V 2 and may 
be strongly reduced by the shock modification (see e.g. Drury et al. l2008h . 

X-ray emissivity due to the thermal bremsstrahlung is given by (Rybicki & Lightman [l979h 

- iSk (ssk) exp (-§) . (C2) 

Here Z is the charge number of plasma ions, g ~ 1 is the Count factor. Assuming the radial density profile 
downstream of the shock p oc r 3< - <7 ~ 1 ^ and the uniform electron temperature we find the flux of the thermal 
X-rays 



2-5 • 10-3^n^| c «| ex p (-f ) . (OS) 

The ga mma- ray flux from the pion decay at energies below the cut-off may be estimated as (e.g. Zi- 
rakashvili [2008h 



ppV ; £ 2 ln(p max /mc) V scm 2?cr pc fcpc ff ^ • 10 3 km s~ V 

(C4) 

Here K„ n = 0.17 is the fraction of the proton energy transmitted to the parent neutral pions, <j pp is the 
total inelastic p — p cross-section, p max is the maximum momentum of accelerated protons and £ cr is the 
ratio of the cosmic ray pressure downstream of the shock P cr to the dynamical pressure PoVf ■ We use the 
value (T pp = 4 • 10~ 26 for this estimate of the flux of TeV gamma-rays. It was assumed that the accelerated 
protons with power-low E~ 2 distribution fill the remnant uniformly. 

The intensity of gamm a-ray s from the IC scattering in the synchrotron losses dominated case may be 
estimated as (Zirakashvili 2008) 



p2 F _ 3£ er K ep R 2 f U rad V f _ erg 2 V f K ep £ er U rad 

IC 8£ B D 2 ln(p max /mc) ° s cm 2 3 • 10 3 km s" 1 10~ 2 & 4 • lO^erg cm" 3 

Here £b is the ratio of the magnetic energy B 2 /8ir to the dynamical pressure PqV 2 , U ra d is the energy density 
of the scattered photons and Rg(deg) is the angular radius of the shock in degrees. 
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The ratio of the gamma-ray energy flux to the ther mal br emsstrahlung energy flux at its maximum at 
Ex = T e is given by equation (see also Katz & Waxman 12008 ) 



E 2 F 



1.76 



6a -3 

CT 2 T l/2 
" ± keV 



Vf 



3- 10 3 km s" 1 



(C6) 



These fluxes are comparable for cosmic-ray modified shocks with £ cr ~ 0.5 and a ~ 6. 



D. Projection effect for a forward and a reverse shock 

The surface brightness J is determined by the integral of emissivity e along the line of sight. For spherically 
symmetric case it can be written as 



oc 

J(r) = 2 J 



r'dr'e(r') 



(Dl) 



For a forward shock the emissivity downstream can be written as e = e$H{R — r) exp((r — R)/l). Here 
H(r) is the step function. For a reverse shock the emissivity has the form e = eoH(r — R) exp((i? — r)/Z). 
For a thin shell of X-ray emission I << R the surface brightness J(r) = y/2Rleob(x) may be written in terms 
of the function b(x) of argument x = (r — R)/l. 

For the forward shock this function (see also Ballet 2006) 



j.) = J 2exp(a;) / dyex.p(y 2 ), x < 



0, x > 



(D2) 



while for the reverse shock 



b(x) — exp(— x) 



2 J dy exp(-y 2 ) , x < 
x/tt, s > 



(D3) 



The function b(x) is shown in Fig. 17. For the forward shock the function b(x) has a width w — 4.6/ at 
the ha lf of m aximum and w = 7.5 1 at th e 1/e of maximum. These numbers are in agreement with results of 
Ballet |2006| and Berezhko & Volk 12004 For the reverse shock this function has a width w = 1.31 at the half 
of maximum and w — 2.6/ at the 1/e of maximum. 
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